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Abstract 

We consider the two-dimensional self-avoiding walk (SAW) in a simply connected 
domain that contains the origin. The SAW starts at the origin and ends somewhere 
on the boundary. The distribution of the endpoint along the boundary is expected to 
differ from the SLE partition function prediction for this distribution because of lattice 
effects that persist in the scaling limit. We give a precise conjecture for how to compute 
this lattice effect correction and support our conjecture with simulations. We also give a 
precise conjecture for the lattice corrections that persist in the scaling limit of the A-SAW 
walk. 
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1 Introduction 



Let D he a bounded, simply connected domain in the complex plane that contains 0. We 
are interested in the self-avoiding walk (SAW) in D starting at the origin and ending on the 
boundary of D. It is defined as follows. We introduce a lattice with spacing 6 > 0, e.g., 67?. A 
self-avoiding walk is a nearest neighbor walk on the lattice with the property that it does not 
visit any site more than once. To be precise, let W(-D, S) denote the set of functions of the form 
w : {0, 1, 2, • ■ ■ , n} — )■ 6Z'^ where n is a positive integer; | \oj{i) — u{i — 1)| | = 6 for i = 1, 2, ■ ■ • , n; 
co{i) 7^ u{j) ioT < i < j < n; uj{0) = 0; E D,j < n; uj{n) ^ D. The integer n is the 

number of steps in the SAW, and from now on we will denote it by |ci;|. 

Since VV(-D, 5) is finite, we can define a probability measure on W(-D, 5) by taking the 
probability of uj to be proportional to where /3 > is a parameter. So 

= VoA^) = ^ (1) 

where the partition function Z{D) = Zs{D) is defined by the requirement that this be a 
probability measure. One can consider this model for all /3 > 0, but it is most interesting for 
one particular value that makes the model critical, (3 = 1/fi, where /z is the connective constant 
which we define next. 

Let Cat be the number of SAW's in the lattice with steps that start at 0. (They are not 
constrained to lie in D.) It is known that this number grows exponentially with A^ in the sense 
that the following limit exists [T3] . 

/i = lim cli^. (2) 

The connective constant fi depends on the lattice. Nienhuis [H] predicted that for the hexagonal 
lattice fi = \/2 + a/2, and this was recently proven by Duminil-Copin and Smirnov For 
the square and triangular lattices there are only numerical estimates of the value of /x. For the 
remainder of this paper we will take /3 = to make the model critical. 

We also consider the analogy of the above definition with the ordinary random walk. The 
natural way to describe the random walk in D starting at and ending on the boundary of 
D is to start a random walk at and run it until it hits the boundary. Let S{D,6) be the 
set of such walks. The probability of a particular such random walk u is 2"'"^' where z is the 
coordination number of the lattice, e.g., z = 4 for the square lattice. We can consider this as a 
random walk starting at the origin stopped when it leaves D, or just as a measure on S{D, 6) 
that assigns measure to each walk of length n. In the SAW the connectivity constant /i 
plays the role of the coordination number z for the random walk. 

In both the SAW and the ordinary random walk we are interested in the scaling limit in 
which the lattice spacing 6-^0. Let us first discuss the case of ordinary random walk which is 
well understood. 
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For the ordinary random walk the scahng hmit is Brownian motion starting at and stopped 
when it hits the boundary of D. The distribution of the endpoint of the Brownian motion on 
the boundary is harmonic measure. The lattice effects associated to the definition of the first 
boundary point of the lattice walk disappear in the scaling limit. The key fact is that if a 
random walk or a Brownian motion gets very close to the boundary, then it will hit it soon. 
Therefore, if we couple a random walk and a Brownian motion on the same probability space 
so the paths are close, then the first time that the random walk hits the boundary will be close 
to the first time that the Brownian motion hits the boundary. See [HI Section 7.7] for a precise 
statement. The argument there works for any simply connected domain even with nonsmooth 
boundaries, and extends to finitely connected domains as well. One does need to assume that 
the boundary is sufficiently large so that when the Brownian motion or random walk gets close, 
then it is very likely to hit it soon. 

If the boundary of our domain D is a piecewise smooth curve, then harmonic measure is 
absolutely continuous with respect to arc length along the boundary [15]. We let hniz) denote 
its density with respect to arc length; this is often called the Poisson kernel (starting at 0). 
If / is a conformal map on D that fixes the origin and such that the boundary of f{D) is 
also piecewise smooth, then the conformal invariance of Brownian motion [12] implies that 
the density for harmonic measure on the boundary of f{D) is related to the density on the 
boundary of D by 

hniz) = \nz)\hf^n)ifiz)). (3) 

If D is simply connected and we take go to he a. conformal map of D onto the unit disc which 
fixes 0, then by symmetry hf(^D-){f{z)) is just l/27r, and so h£){z) = \g'^{z)\/2n. We emphasize 
that conformal invariance of harmonic measure implies that the exponent of in ([3]) equals 

one. 

We now consider the SAW in D from the origin to the boundary of D. Here we will assume 
that dD is piecewise smooth. In this case we have the following conjecture. 

• As 5 — )• 0, the measures Pd,<5 converge to a probability measure Pd on simple paths from 
the origin to dD. 

• The measure Pn can be written as 

/ pD{z)PD{0,z)\dz\, (4) 

JdD 

where Pd{z) is the density of a probability measure on dD and P_d(0, z) is the probability 
measure associated to radial SLE^/s from to z in 

• There exists a periodic function, l{6), such that 

pn{z) = mz,D))pn{z) (5) 

where 9{z, D) is the angle of the tangent to dD at z and Pd{z) is a multiple of the SLEg/^ 
partition function. The function l{6) and its period depend on the lattice. For example, 
on the square lattice the period is ti/2. 
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• The density Pd{z) transforms under conformal maps by 

pn{z) = c\f{z)f'pf^n)U{z)) (6) 

The constant c is determined by the constraint that this be a probabihty density. (It 
depends on D and /.) In particular, if D is simply connected, and qd '■ D ^ with 
^^(0) = 0, then 

pn{z) = c\Uz)\'/' (7) 
(D denotes the unit disc centered at the origin.) 

In the case of 1? where the boundary of D is composed of horizontal and vertical line 
segments, this conjecture was made by Lawler, Schramm and Werner [lOj. Simulations on an 
infinite horizontal strip [3] give strong support to the conjecture. For such boundaries, l{6{z, 0)) 
is constant. The conjecture was reiterated in [9] where it was also conjectured for other domains 
"...after taking care of the local lattice effects." Our conjecture makes precise the nature of the 
lattice effects in terms of the lattice correction l{6{z, 6)). The function l{6) will depend on just 
how we define "ending on the boundary of and on the lattice type. In general, we use p to 
denote densities that do not include the lattice effects that persist in the scaling limit, e.g., (|6]), 
and we use p to denote densities that do include the lattice effects, e.g., ([5]). 

Unfortunately, the Monte Carlo algorithms for simulating the above ensemble are local 
algorithms and so are not very efficient [13]. We have not attempted to test the conjecture for 
this ensemble by simulation. Instead we introduce another ensemble that we can study with the 
pivot algorithm, a fast global Monte Carlo algorithm. Instead of stopping a walk at a boundary 
point, one chooses an infinite length walk conditioned on the event that it crosses the boundary 
only once. We will refer to this ensemble as the "cut-curve" ensemble since the boundary of 
the domain cuts the SAW into two SAW's, one contained in D and one in the complement of 
D. The above ensemble and the cut-curve ensemble are the same for the infinite strip studied 
in [3], but for most domains this is not the case. 

The scaling limits of the SAW and the loop-erased walk, which is obtained by erasing loops 
from the ordinary random walk, are two cases of the Schramm-Loewner evolution (SLE). The 
discrete models can be considered as special cases of the A-SAW. We review the A-SAW in 
the next section, and we extend our conjecture to this case. (There is no precise conjecture 
in the literature on the nature of the lattice correction, and we think it is worthwhile to write 
it down.) After that we consider the conjectured scaling limits of the two ensembles — walks 
stopped upon reaching the boundary and walks conditioned to hit the boundary only once. In 
only the k = 8/3 (SAW) case do we expect an equivalence of these ensembles, and therefore 
the tests we do here would not work for other values of /t. 

In section two we return to the SAW and give explicit conjectures for the lattice correction 
function l{9) for the two lattice ensembles. In section three, we discuss simulations for the 
cut-curve ensemble including numerical calculation of the lattice correction function. In the 
final section we summarize our results and discuss the lattice correction function l{9) for other 
interpretations of the SAW ending on the boundary of the domain. 
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1.1 A-SAW 



The conjectures for the SAW and the results for the ordinary walk (considered in terms of the 
loop-erasure of the paths) are particular cases of conjectures for a model called the A-SAW 
introduced by Kozdron and Lawler [6]. Even though we are only testing the SAW conjecture, 
we will give the conjectures for the general model. There are two versions, chordal (boundary- 
to-boundary) and radial (boundary-to- interior); we will restrict our discussion here to the radial 
case. (There is also an interior-to-interior case, but then there is no boundary lattice correction 
so it is not relevant for this paper.) As above, we assume that D is a domain with piecewise 
smooth boundary. For convenience we use 1? for our lattice, but the definition can be extended 
to other lattices. The parameter A can be considered a free parameter, but we will set A = — c/2 
where c < 1 denotes central charge. (We will not define central charge in this paper and can 
just take it as a parameter.) The A-SAW is a model that is conjectured to have a scaling limit 
of SLE^ where 

13 -c - v/(13 - c)2 - 144 , , 
K = e (0,4]. 

The (rooted) random walk loop measure is the measure on ordinary random walk loops 
that assigns weight (1/4)" to each loop of > steps. A loop is a path which begins and 
ends at the same point. The measure can also be considered as a measure on unrooted loops 
by forgetting the root. 

For each lattice spacing S and each SAW co as above, we let itid^uj) denote the total measure 
of the set of loops that lie entirely in D and share at least one site with u. If /3 > 0, the A-SAW 
gives each SAW u as above weight 



g(t^) = q5,D,cA^) = 1^ exp |--mD(w) 
The partition function is 

Z{D;6) = Z{D;6,c,^) = J2<li^)^ 

where the sum is over all SAWs on the lattice that start at the origin and end at the 
boundary. (As noted before, there are several definitions for "ending at the boundary" . In the 
discussion here, we fix one such definition and the lattice-dependent quantities depend on the 
choice.) It is conjectured that for each c, there is a critical value /? = /3c such that Z{D; S) 
follows a power law in 5 as 5 — )■ 0. We assume this conjecture and fix /3 at the critical value. 

If c = 0, then this is the usual SAW model. If c = —2, then this is the loop-erased random 
walk (LERW) which is obtained by taking the ordinary random walk as above and erasing 
loops chronologically from the paths. The partition function Z{D;6) for the loop-erased walk 
is exactly the same as that for the usual random walk. (See [8], Chapter 9] for a discussion 
of this.) We state our conjectures in terms of the boundary and interior scaling exponent for 
SLE,: 

6 - K ,^ , ~ 6 (k - 2) 
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• There is a lattice correction function l{6) that is continuous, strictly positive and periodic. 

• If is from to dD, let l{u) denote l{0{z, D)) where z is the first point on dD hit by a 
bond of u and 9{z, D) is the angle of the tangent to dD at z. Define the lattice-corrected 
weight by 

= WY 

• As (5 — )■ 0, the measure ^^^^^^ q on paths converges to a nontrivial finite measure vd on 
simple paths from to dD. It can be written as 

T^D= I Pd{z) vti^.z) \dz\, 

JdD 

where po is a positive function and z/*(0,2) is a probability measure on simple paths 
starting at and leaving D at z. 

• If (yf is a conformal transformation with g{0) = that is smooth on dD, then 

PD{z) = \9'{z)\'\g'm'p,^o){z). (8) 

• The probability measures ^'*(0, z), considered as measures on curves modulo reparametriza- 
tion, are conformally invariant. More specifically, if 7 : [0,to) — )■ -D is a curve with 
7(to-) e dD, let 

a{t) = j\g'{^{s))\Us, d=l + ^, 

and define g o •jit) by 

goj{cr{t))=g{j{t)), 0<t<a{to). 
li g o /i*(0, z) denotes the induced measure on curves on g{D), then 

gou*{0,z) = z/^%)(0,2;). 

• If D is simply connected, then i/*(0,2;) is the reversal of radial SLE^^ from 2; to in D 
with the natural parametrization. (This should also be true for multiply connected D 
under the appropriate definition of SLE^ proposed in [7].) 

These conjectures are a long way from being proved. Indeed, a special case is the SAW 
model which is a notoriously difficult problem! However, we state them here to see that the 
precise conjecture requires discussing a boundary lattice correction; our conjecture is that the 
correction only depends on the angle of the boundary. We make a number of comments. 



6 



• The density Pd{z) is sometimes called the partition function for radial SLE,^. It is defined 
up to an arbitrary constant. 

• The loop-erased walk (c = —2) is particularly nice because it is closely related to the 
ordinary random walk. The partition function is the Poisson kernel (density of harmonic 
measure) even for non-simply connected domains. 

• For other values of c, if D is simply connected, the partition function is the Poisson kernel 
raised to the 6th power. However, this is not true for multiply connected domains. 

• The loop terms m£){ijj) can be considered as having three parts. The very small loops 
that occur away from the boundary contribute a microscopic (lattice dependent) part 
that affects the critical value of /3. The large loops contribute a macroscopic term that 
is seen in the scaling limit; this is the Brownian loop measure as defined in [TT]. Finally, 
there are the small loops that occur near the boundary. They give both a macroscopic 
effect seen in the exponent b and a microscopic effect in the function /. The boundary 
effect is measured both in the number of walks that stay one one side of a line and in the 
measure of loops that go on the walks. We only see the first effect in the SAW (c = 0) 
case. 

• If K = 8/3, then h = 5/8, b = 5/48. There is a |5''(0)|'' factor in ([8]) that does not seem to 
appear in (Q. However, it is implicitly there in the normalization to make the measure a 
probability measure. 

1.2 Cut-curve configurations 

We will be testing a cut-curve ensemble for SAW's. The case c = 0, which is what we use in this 
paper, is special for such configurations and agrees with the bridge decomposition of restriction 
measures p], but for the sake of completeness let us discuss the general case for c < 1. Suppose 
D is a bounded, simply connected domain containing the origin whose boundary dD is a smooth 
Jordan curve. Let D* be the unbounded component of C \ dD. We consider two measures on 
simple curves from to infinity that intersect dD only once: 

1. Take "whole-plane" SLE^ and condition on the event that the curve intersects dD only 
once. This is conditioning on an event of probability zero, so one must define this in 
terms of a limit. 

2. Take independent copies of radial SLE^ in D and in D* and condition them to hit the 
same point in dD. Then concatenate the paths. 

In order to see the difference, let us describe the lattice models that we expect to converge to 
these measures. Since it is hard to talk about infinite walks, we will choose a point z* G D* 
with large absolute value. We take the scaling limit for fixed z* and then let z* go to infinity. If 
(5 is a scaling factor, we abuse notation slightly and write z* for the point in closest to z*. 
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We consider two measures on paths. For each S, we consider the set W(5, D, z*) of SAW's u 
on (5Z^ with uj{0) = 0,uj{\uj\) = z*, \u){j)\ < \z*\'^ for all j such that there exists only one bond 
that intersect dD. We write z^j, w^j for the vertices in this bond with z^^ E D,Wi^ E C \ D, and 
we write u = © [z^^, w^j] © co* where 

= [u{0), . . . ,Zu;], u* = [w^,. . . 

Let (3 = (3c he the critical value. Then we consider the following measures on W((5, D, z*): 

1. Each oj gets weight 

qi{uj) = Z?!"^! exp |-|m(w; , 

where m{u; z*) denotes the measure of loops staying in the ball of radius \z*\'^ that 
intersect u. 

2. Each OJ gets weight 

q^iu) = /3H exp {-£[mz,K) + mD^u*)]} 

We can state our conjectures as follows. Let 9{u) be the angle of the tangent to dD at the 
point where u hits dD. 

1. There exists a lattice correction function li such that we can take the scaling limit of the 
measure 

^i(^) = 77T- 

If we then take z* — >■ oo, we get whole plane SLE,^ conditioned to hit dD only once. 

2. There exists a lattice correction function I2 such that the scaling limit of the measure 

exists. If we then take z* to infinity, we get the measure given by 

Pd{z) p*i){z) z/^(0, z) © 00) \dz\. 

IdD L ^ 

Here p*{z) is the density as in (j8]) for radial SLE^^ in D* centered at infinity. 

For c = (restriction measures), the two limits agree and this is what we use for SAW in 
this paper. For other values of c we get different measures. For example if c = —2 (loop-erased 
walk), the first measure corresponds to loop-erased walk conditioned to hit dD only once and 
the second measure corresponds to the loop-erasure of an ordinary walk conditioned to hit dD 
only once. 
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2 Lattice effects 



The constraint that the SAW stays in D has both a macroscopic and a microscopic effect on the 
boundary density. The macroscopic effect is captured by the conjecture (EI). The microscopic 
effect comes from the behavior near the endpoint of the walk on the boundary of D. Consider 
a SAW that ends at z G dD and consider the tangent hne to the boundary at z. The constraint 
that the SAW stays in D imphes that near z the SAW must stay on one side of this hne. Loosely 
speaking, the number of SAW's that end at z and stay on one side of this line depends on the 
orientation of the line with respect to the lattice. The result is a factor l{6) that depends on 
the angle of the tangent line with respect to the lattice, and so we obtain our conjecture ([5]). 

The lattice correction function l{6) depends on the type of lattice and on how we interpret 
"ending on the boundary of D." We will first discuss the interpretation we introduced at 
the start of the introduction. We consider all SAW's u with ci;(0) = 0, u^i) E D ioi i = 
0, 1, 2, ■ ■ ■ , 1^1 — 1 and co{\oj\) ^ D. (|w| denotes the number of steps in co.) So the last bond of 
the SAW intersects the boundary of D, and this is the only bond in the SAW that intersects 
the boundary. 

We can compute the lattice correction function l{9) as follows. We give the details for the 
square lattice. Other lattices, e.g., the triangular or hexagonal, will require some modifications. 
Consider a bond that intersects the boundary dD. Let z be the endpoint of the bond that is 
in D. Let w be the point where the bond intersects the boundary (typically not a lattice site). 
Consider the tangent line to the boundary at w. We need to count the number of SAW's of a 
fixed length that end at z and do not intersect this tangent line. This quantity will depend 
on the angle of the tangent line with respect to the lattice. It will also depend on the distance 
|w — 2;| and on whether the bond is horizontal or vertical. However, these two factors are the 
only dependence on the bond. 

Motivated by the above, we consider the vertical bond between (0, 0) and (0, 1). We fix an 
/ in the interval [0, 1] and an angle 6. The parameter / plays the role of the distance |w — 
Let L be the line with polar angle 6 passing through the point (0, /). We consider SAW's with 
N steps ending at the origin which by reversal can be considered as beginning at the origin. 
Let ctv be the number of such walks, and let a]\f{l,6) be the number of such walks that do not 
intersect the line. So a^il, 0)/cn is the probability that the step SAW does not intersect the 
line. We expect that there exists a function 6) such that 

Cn 

with p = 25/64 [lOj. The actual exponent is not important here. The key is that if I', 6' are 
two different values, then 



We define 



A^ ^ oo. 



K(/,^)= hm ^^AT^ (9) 
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(We do not know how to prove this hmit exists.) The superscript f on p indicates that we 
took the bond crossing the boundary to be vertical We use the subscript 1 to distinguish 
the quantities related to the lattice correction function for this particular ensemble from the 
ensemble that we will consider next. We let Pi denote the analogous quantity for a horizontal 
bond. The symmetry of the square lattice implies that Pi{l,9) = p\{l,9 + vr/2). (Note that 
Pi{l, 6) and p\{l, 6) have period tt in 6.) 

Now consider what happens as we move along the boundary. The angle 9 of the tangent will 
vary smoothly. The distance / will not. As long as tan(0) is not rational, / will be distributed 
uniformly between and 1. So in the scaling limit, averaging over an infinitesimal section of 
the boundary will be equivalent to averaging / over [0, 1]. So we define 

pl[e) = [ pUi,e)di, x = h,v 

Jo 

The function Pi{6) captures the microscopic lattice effect caused by the constraint that 
the SAW must stay on one side of a line as it approaches the boundary. There is another 
lattice effect that comes from the density of bonds that cross the boundary. Define b^{d) to 
be the density of vertical bonds along a line with polar angle 6, i.e., the average number of 
vertical bonds that intersect the line per unit length. We have b^{9) = \ cos(6')|. The density of 
horizontal bonds b'^{6) is b"{6 — 7r/2) = \ sin(^)|. The lattice correction function is then 

k{e) = b''{e)piie) + b\e)p',{e) (lo) 

and the boundary density for this first interpretation of ending on the boundary is 

pnA^) = c\\gU^)\'/%ie{z,D)) (11) 

The Monte Carlo algorithms for simulating the above ensemble are local algorithms and so 
are not very efficient [13]. We have not attempted to test the conjecture for this ensemble by 
simulation. Instead we use the cut-curve ensemble introduced earlier. We consider the infinite 
length SAW in the full plane starting at 0. We condition on the event that the SAW crosses the 
boundary of D exactly once. Consider a bond crossing the boundary, and let z be the endpoint 
inside D and w the endpoint outside of Z^. If we condition on the event that the SAW contains 
this particular bond and this is the only bond in the SAW that crosses the boundary, then we 
have a SAW in D from to z and a SAW in the exterior of D from w to oo. So the conjecture 
for the boundary density will be a product of two functions. The interior SAW from to z 
gives a factor of Pd{z) with Pd{z) given by ([6]). The exterior SAW from w to oo gives a factor 
of p*£,{z) given by 

pUz) = c\h'^{z)\'/' 

where hj:,{z) is the conformal map of D* onto the unit disc D with hoioo) = 0. Our conjecture 
for the density of the point where the SAW crosses the curve is 

PD,2{^)=PD{z)pU^)l2{0iz,D)), (12) 
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where hi^) is the lattice correction function for the cut-curve ensemble. (The subscript 2 
indicates that the quantities are for the cut-curve ensemble.) 

We can compute hi^) as follows. Again, we restrict our attention to the square lattice. For 
/ G [0, 1] and an angle 9, let L be the line with polar angle 9 passing through (0, /). We consider 
SAW's with N = 2n + l steps such that the middle bond is the bond between (0, 0) and (0, 1). 
Let d]sr be the number of such walks, and let bi\f{l,9) be the number of such walks that only 
intersect the line once. (Of course it must be the middle bond that has the intersection). (By 
translating our SAW's so that they start at 0, we see that = ciq/2.) Note that biy{l,9)/dN 
is the probability the N step SAW has no intersection with the line other than the middle 
bond. We can think of generating this step SAW by first generating two n step SAW's which 
are independent and attach to (0, 0) and (0, 1) and then keeping them only if they mutually 
self avoid. The probability they are mutually self- avoiding goes as A^^"''' with 7 = 43/32. The 
probability that both of the n step SAW's do not intersect the line goes as N~'^p with p = 25/64 
[lOj. Note that if both of them do not intersect the line, then they are mutually self-avoiding. 
Hence we expect that there exists a function P2{1, 9) such that 

dn 

We define 

pl{l,9)= hm Ml^iv2p-7+i. 

' N^oo dM 

Then we define 

pm= l\i{i.o)di 

Jo 

As in our derivation of the lattice correction function for the first ensemble, the integral over I 
comes from averaging over bonds crossing a small section of the boundary. We define p^{l,9) 
and P2(^) i^ analogous way. As before we let h^{9) and h"{9) denote the densities of 
horizontal and vertical bonds along a line with polar angle 9. The lattice correction function is 
then 

l,{9) = h\9)pl{9) + h\9)p^,{9) (13) 

In the above discussion we have considered SAW's starting at an interior point in the domain. 
The same discussion applies to an ensemble of SAW's that start at a prescribed boundary point 
of D. Th conjecture for the boundary density again transforms according to (jH]). With the 
starting point on the boundary this density is not normalizable. We must restrict the endpoint 
of the SAW to a subset of the boundary that is bounded away from the starting point to get a 
normalizable density. A useful reference domain in this case is the upper half plane with the 
starting point at 0. The unnormalized density for the harmonic measure is and so for the 
SAW it is x-^l^. 
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3 Simulations 



In this section we study the cut-curve ensemble by Monte Carlo simulations. There are two 
types of simulations. We compute the lattice correction function hiO) by simulation, and we do 
simulations of the SAW in two different geometries to test the conjecture (fT2|) . We first discuss 
the computation of hiO). 

Recall that for odd N, dN is the number of SAW's with N steps such that the middle bond 
is the bond between (0,0) and (0, 1). For / G [0, 1] and an angle 9, bp^{l,9) be the number of 
such SAW's whose only intersection with the line through the point (0, 1) at angle 6 is in this 
middle bond. The ratio b^{l,9)/diy is a probability and so may be computed as follows. We 
use the pivot algorithm to generate SAW's with steps that start at the origin and such that 
the middle bond is vertical. We then pick a point on this bond uniformly at random and take 
the line with angle 6 to go through this point. We test if the only intersection of the SAW 
with the line is through the middle bond. We find the fraction of the samples that satisfy this 
condition and multiply it by A^2p-7+i^ fjj^g result is an estimate of p^i^)- Note that we have 
included the integral over / from to 1 in the simulation by randomly choosing / uniformly 
from [0, 1] for each sample. 

We did this simulation for values of N ranging from 101 to 5001. For the smaller values 
of A^, one can clearly see finite A^ effects. As is always the case with simulations of the SAW, 
the time required grows with A^. However, in this simulation this is exacerbated by the fact 
that the probability the SAW only intersects the line once goes to 0, and we must multiply 
the probability we compute in the simulation by A^2p-7+i Even with 2 billion samples, our 
simulations for A^ = 2001 and 5001 have significant statistical errors. The results for A^ = 501 
and A^ = 1001 differ by at most 0.05%, and we use A^ = 1001 in this paper. We generated 5 
billion samples for this case which took about 67 cpu-days on a rather old cluster with 2.4 GHz 
cpu's. Figure [1] shows the function p^i^)- figures we give the angle 9 in degrees. 

We now turn to the simulations to test conjecture (fT2|) . We use the pivot algorithm to 
generate walks in the full plane with a constant number of steps A^. We take A^ = 1, 000, 000. 
The lattice spacing and A^ are such that the size of the SAW is large compared to the domain 
D so that the SAW is effectively infinite. We condition on the event that the SAW intersects 
the boundary of D exactly once. The probability of this event goes so zero in the scaling limit, 
so we must generate very large numbers of SAW's to get good statistics. We use Clisby's fast 
implementation of the pivot algorithm using binary trees p]. 

In our first test we take the domain for the cut-curve ensemble to be a disc centered at 
where the SAW starts. We take the lattice spacing to be and take the radius of the disc to 
be i? = 0.2. (With R = 0.3 the effect of the finite length of the SAW begins to be noticeable. 
At R = 0.4 it is quite noticeable.) In the simulation we sample the Markov chain every 1000 
iterations and generate a total of approximately 47 million samples. Just over 10% of these 
samples satisfy the condition that the SAW only intersects the boundary of the circle once, and 
we have approximately 4.9 million samples of the boundary density. 

In this geometry both Pd{z) and p*£,{z) are constant, so the prediction for the boundary 
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Figure 1: The function p^iO)- 

density without lattice effects is just the uniform density. The angle 6{z, D) of the tangent line 
at z is equal to the polar angle 6 oi z mod 90 degrees. So if we think of the boundary density 
as a function of the polar angle 6, then our conjecture ffT^ is that the boundary density is 
proportional to hiO). Figure [2] shows the function hiO) and the boundary density we find in 
the simulation of the cut-curve ensemble. Both functions have a period of 90 degrees. We plot 
the boundary density as a function of 9 mod 90. Both curves are normalized so that the total 
area under the curves is one. 

Figure [2] compares densities. Actually, the function plotted for the simulation of the cut- 
curve ensemble is a histogram. So the points plotted correspond to the average value of the 
density over a small interval. The simulations do not compute densities directly. Finding the 
density requires taking a numerical derivative, i.e., computing a histogram. We can avoid this 
extra source of numerical error by working with cumulative distribution functions (cdf 's) rather 
than densities. 

In figure [3] we study the cdf for the cut-curve ensemble using the unit disc centered at 0. We 
plot two curves. One is the cdf we find in the simulation of the cut-curve ensemble minus the 
cdf for the uniform density, i.e., the density given by ([7]). (In this figure we have again taken 
advantage of the periodicity of the underlying density functions.) This difference is small with 
the maximum being slightly less than 2%, but it is clearly not zero. In the second curve we show 



13 



0.014 




0.01 - 



0.009 - 



0.008 I 1 1 — lJ 1 l_iJ 1 Lj_J 1 li — I 1 — il 1 l__iJ 1 1 Li — I 1 1 

10 20 30 40 50 60 70 80 90 

Theta 

Figure 2: The function l2{9) (red curve) and the density for the full plane SAW 
conditioned to hit a circle exactly once (blue histogram). Note that the range of the 
vertical axis does not start at 0. 

the cdf for the simulation of the cut-curve ensemble minus the cdf corresponding to the density 
given by f[T^ . i.e., corresponding to hid)- The difference is on the order of 0.02%. The error 
bars in the figure are two standard deviations for the statistical errors, i.e., the error that comes 
from not running the Monte Carlo simulation forever. There are also errors in the simulation 
from two other sources - the finite length of the SAW and the nonzero lattice spacing. We have 
studied the error from the finite length of the SAW by simulating the ensemble with several 
values of the radius of the disc. With R = 0.2 we believe that the error from the finite length 
of the SAW is much smaller than the statistical errors. The nonzero lattice spacing means that 
all our random variables are at a small scale discrete random variables. This is reflected in the 
slightly jagged nature of the second curve. The error from the nonzero lattice spacing appears 
to be comparable in size to the statistical error. Thus the difference between the cdf from the 
simulation and the cdf given by f|T2|) is zero within the errors in our simulation. Figure [3] gives 
evidence that there are indeed lattice effects that must be taken into account in the boundary 
density and that our conjecture correctly accounts for these lattice effects. 

For our second test of conjecture (fT2|) . we consider the SAW in the upper half plane, starting 
at the origin. We take the cut-curve to be a semi-circle centered at the origin. Again, we take 
the lattice spacing to be N~'^ and the radius of the semi-circle disc to be i? = 0.2. We sample the 
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Figure 3: The simulation of the cdf for the boundary density for the cut-curve en- 
semble for the first geometry minus the theoretical prediction. The larger curve does 
not include the lattice correction function; the smaller curve does. 

Markov chain every 1000 iterations and generate a total of approximately 27 million samples. 
Approximately 13% of these samples satisfy the condition that the SAW only intersects the 
boundary of the circle once, and we have approximately 3.6 million samples of the boundary 
density for this geometry. 

The ensemble consists of all SAW's in the upper half plane which start at and only cross 
the semicircle once. In this geometry the arc length along the semicircle equals the polar angle 
9. So we will express densities as functions of 9. A simple computation using the conformal 
map z + 1/z shows that the interior density Pd{9) is [sin(^)]^/*. The exterior density p*d{9) 
is exactly the same. (This is just a consequence of the symmetry of our geometry under the 
inversion z — )■ —1/z.) So our conjecture for the density along the cut-curve is proportional to 

[sm{9)f^ l{9) (14) 

The comparison of our simulation of the cut-curve ensemble cdf for the SAW and our 
conjecture is shown in figure HI Again we plot two curves. For one curve we find the cdf 
corresponding to just the density function [sin(6')]^/^. This would be the conjectured cdf if 
there were no lattice effects. We plot the cdf for the simulation of the cut-curve ensemble 
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Figure 4: The simulation of the cdf for the boundary density for the cut-curve en- 
semble for the second geometry minus the theoretical prediction. The larger curve 
does not include the lattice correction function; the smaller curve does. 

minus the cdf corresponding to just [sin(6')]^/^. This difference is small (the max is on the order 
of 1.5%), but is clearly not zero. For the other curve we compute the cdf corresponding to our 
conjecture f|T4|) with the lattice effect and subtract this function from the simulation of the cdf 
for the cut-curve ensemble. The difference is on the order of 0.05% which is zero within the 
errors in our simulation. This figure gives further evidence that there are indeed lattice effects 
that must be taken into account in the boundary density and that our conjecture correctly 
accounts for these lattice effects. 

4 Conclusions and future work 

We have considered the ensemble of SAW's in a simply connected domain containing the origin 
which start at the origin and end on the boundary. It has been noted before that the conjecture 
for this boundary density will have lattice effects that persist in the scaling limit. We have 
conjectured that this lattice effect is given by multiplying the density by a function l{9{z,D)) 
where 6{z, D) is the angle of the tangent line to the boundary of D at the point z G dD. The 
lattice correction function l{6) depends on the lattice and on how we interpret "ending on the 
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boundary of DP We have shown how to compute the lattice correction function l{9) for two 
particular interpretations. Our focus has been on the distribution of the endpoint of the SAW 
on the boundary, but we should remark that in light of (jlj) the lattice effects in this boundary 
density will produce lattice effects in the probability measure Pd- We have also extended this 
conjecture to the A-SAW. 

As we have noted before, there is no efficient way to simulate the natural interpretations of 
the ensemble of SAW's in a domain which start at an interior point and end on the boundary. We 
have circumvented this difficulty by introducing the cut-curve ensemble which can be thought 
of as an ensemble of two SAW's, one from the interior point to the boundary and the other from 
that boundary point to oo. There is another ensemble that is amenable to efficient simulation 
which is studied in [5j. Given a domain D containing the origin the ensemble is defined as 
follows. We assume the domain has the property that a ray from the origin only intersects 
the boundary of the domain in one point. For a SAW u we let X{u)) > be such that the 
endpoint of u is on the boundary of X{uj)D. In general the SAW need not be inside the dilated 
domain \{u)D. Our ensemble consists of all SAW's uj of any length such that uj is contained 
in \{uj)D. (One must introduce cutoffs to make this a finite measure.) In [5] we show how one 
can simulate this ensemble using the ensemble of SAW's of a fixed length. 

Finally, it is natural to ask if there is an interpretation of "ending on the boundary of D" for 
which there are no lattice effects in the scaling limit. We speculate that the following ensemble 
has this property. As before, let 5 be the lattice spacing. Let e > and consider all SAW's that 
start at the origin, stay inside D and end within a distance e of the boundary of D. We let 5 
go to first and then let e go to zero. We conjecture that l{9) is constant for this ensemble. 
Unfortunately, the double limit involved in this ensemble makes it difficult to simulate this 
ensemble. 
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